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^ : Abstract 

' We propose a collective opinion formation model with a so-called confirmation bias. The confir- 

'. 

^~>' mation bias is a psychological effect with which, in the context of opinion formation, an individual 

; 

in favor of an opinion is prone to misperceive new incoming information as supporting the current 
belief of the individual. Our model modifies a Bayesian decision-making model for single Individ- 

> 

0> ' uals (Rabin and Schrag, Q. J. Econ. 114, 37 (1999)) to the case of a well-mixed population of 

. interacting individuals in the absence of the external input. We numerically simulate the model to 

^SJ \ show that all the agents eventually agree on one of the two opinions only when the confirmation 

o 

en ' bias is weak. Otherwise, the stochastic population dynamics ends up creating a disagreement con- 

figuration (also called polarization), particularly for large system sizes. A strong confirmation bias 



^ , allows various final disagreement configurations with different fractions of the individuals in favor 

of the opposite opinions. 
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I. INTRODUCTION 



There are various models of collective opinion formation in which agents modify their 
opinions according to interaction with other agents . Opinion formation is a dynamic 
process: for example, interaction between agents makes their opinions approach each other. 
An important problem in opinion dynamics is to examine assessing when an agreement (i.e., 
consensus) among all the agents occurs. Complete agreement is rarely observed in the real 
world jsl, 0] . However, it is an established fact that opinion dynamics under the voter model, 
a classical opinion model in statistical physics and probability theory, inevitably reaches 
agreement in finite populations Isl-Ql. The majority rule model has a similar feature 



id, lUj- Partly motivated by this discrepancy, various extensions of voter and majority 



rule models and different models of collective opinion formation have been proposed to 
account for the disagreement in finite populations. Examples include language competition 



models 
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2l-ll9|. voter-like models on adaptive networks 2o|-23|, voter model under partisan 
bias (the assumption that agents naturally prefer one opinion) 2J, |25j, and variations of 
Axelrod's cultural dynamics (see |l| for references). Theoretical models have also been 
proposed in social sciences to explain disagreement in the context of polarization. For 
example, prior beliefs or initially received signals can cause disagreement between agents. 



even if they receive the same public signals from then on [26- 

Although there is a plethora of studies addressing the problem of agreement and dis- 
agreement in opinion dynamics, we propose a new model incorporating two factors that 
are relevant to human behavior: Bayesian belief updating and confirmation bias. Bayesian 
belief updating is commonly used in studies of the decision- making of agents receiving uncer- 
tain information 
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3l| . The confirmation bias is a psychological bias inherent in humans. 



in which an agent inclined towards an opinion tends to misperceive incoming signals as 
supporting the agent's belief 
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33|. 



We are not the first to study opinion formation under the Bayesian updating and confir- 
mation bias. In the framework of single agent opinion formation, Rabin and Schrag showed 
that the confirmation bias triggers overconfidence and can cause the individual to hold in- 
correct beliefs, even if it receives a series of external signals suggesting the true state of the 
world 3J]. Orlean studied the Bayesian dynamics of agents subjected to the confirmation 



bias, interacting through the mean field 



35| . The model yields agreement or disagreement 



depending on the parameter values. 



34| . we propose a model of coUec- 



In this study, motivated by the Rabin-Schlag model 
tive opinion formation with a confirmation bias. We model direct peer-to-peer interactions 
between agents (not through the mean field) and their effects on the Bayesian updating 
of each agent. To study the pure effects of interactions among agents, we do not assume 



that agents receive signals from the environment as in previous studies [30| . We numerically 
simulate the model to reveal the conditions under which the population of agents agree and 
disagree, depending on the values of parameters such as the strength of the confirmation 
bias, fidelity of the signal, and the system size. 



II. MODEL 



Our model modifies the Bayesian decision-making model proposed by Rabin and Schrag 



34j | in two main ways. First, we consider a well- mixed population of Bayesian agents that 
interact with each other; Rabin and Schrag focused on the case of the single agent. Second, 
agents do not receive external signals from the environment in our model. In the Rabin- 
Schrag model, such an external signal, which represents the "correct" answer in the binary 
choice situation (i.e., the true state of nature), is assumed. By making the two changes, 
we concentrate on collective opinion formation by Bayesian agents, whereby there are two 
possible alternative opinions of equal attractiveness. 

We label the agents 1, . . . , and denote the opinion of agent i {i = 1, . . . , A^) by 
Xi G {A,B}, where A and B are the alternative opinions. We assume that agents are not 
perfectly confident in their opinions. To model this factor, we adopt the Bayesian formalism 



used by Rabin and Schrag [3J|. We denote by Pr(xi = A) the strength of the belief, 
(hereafter, simply the belief), with which agent i believes in opinion A. A parallel definition 
is applied to Pr(xj = B). It should be noted that Pr(xj = A) > 0, Pr(xj = B) > 0, and 
Pr(xj = A) + Pr(xj = B) = 1. If Pr(xj = A) = 1/2, agent i is indifferent to either opinion. 

We update the agent's belief as follows. In each time step, we select an agent i to be 
updated with equal probability 1/A^. Agent i refers to agent j's opinion for updating i's 
belief Pr(xj = A), where j (7^ i) is selected with equal probability 1/(A^ — 1) from the 
population. Agent j imparts a signal s G {a, b}, where a and b correspond to j's opinions A 
and B, respectively. We assume that the probabilities that agent j imparts s = a and s = b 



are given by 



Pr(s = a) = Pr(s = a\xj = A) Fi{xj = A) + Pr(s = a\xj = B) Fr{xj = B) 
= e FT{xj =A) + {l-e) FT{xj = B) 

= 1-6 + {26-1) Fiixj = A) (1) 



and 



Pr(s = b) = Pr(s = b\xj = A) Pr(xj = A) + Pr(s = b\xj = B) Ft^xj = B) 
= (1 - ^) Fr{xj = A)+6 FT{xj = B) 

= 6 - {26 -l)Fi{xj = A), (2) 



respectively, where 



6 = Fi{s = a\xj = A) = Pr(s = b\xj = B) (3) 



represents the rehabihty of the signal and 1/2 < 6 < 1. If j is confident in its own opinion 
and the transformation from j's belief (i.e., Pr(xj = A)) to j's output signal (i.e., a or 
b) is reliable, signals a and b are likely to indicate opinions A and B, respectively. In 
the limit 6^1, Pr(s = a) ^ Ft^xj = A) and Pr(s = b) ^ Fi{xj = B). If ^ = 1/2, 
Pr(s = a) = Ft{s = b) = 1/2 such that s does not convey any information about j's belief. 
We implicitly assume that all the agents share the same value of 6 and that they know this 
fact when performing the Bayesian update, as described below. 

When agent j imparts signal s G {a,b}, agent i is assumed to perceive a subject signal 
cr G {a, (3}, where a and (3 correspond to A and B, respectively. The flow of the signal 
conversion is depicted in Fig. [H If agent i is not subject to the confirmation bias, a and f3 
are equal to a and b, respectively. Otherwise, agent i may misinterpret the signal imparted 
by agent j, depending on the prior exposure of agent i to other signals. Following Rabin 
and Schrag [Mi], we define 

Pr((T = a\s = a, FT{xi = A) > 1/2) = Fi{a = (3\s = b, Fi{xi = A) < 1/2) = 1 (4) 

and 

Pr(cT = a\s = b, Pr(xi = A) > 1/2) = Pr(cT = = a, Pr(xi = A) < 1/2) = q, (5) 
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where q {0 < q < 1) parameterizes the strength of the confirmation bias. Equation ([5]) states 
that an agent preferring opinion A misinterprets an arriving b signal as A (i.e., a = a) with 
probabihty q. If g = 0, the confirmation bias is absent, and s = a and s = b are always 
converted to cr = a and a = P, respectively. If g = 1, the agent perceives the signal that is 
consistent with its current preference (i.e., a if Pr(xj = A) > 1/2 and /3 if Pr(xj = A) < 1/2), 
irrespective of the signal imparted by agent j (i.e., a or 6). The other conditional probabilities 
can be readily derived from Eqs. (jl]) and ([5]). For example, Eq. (jlj) implies 

Pr(a = I3\s = a, Fr{xi = A) > 1/2) = 1 - Pr(cT = a\s = a, Pr(xi = A) > 1/2) = (6) 

and Eq. ([5]) implies 

Pr(cT = a\s = a, Pr(x, = A) < 1/2) = 1 - Fi{a = l3\s = a, Fi{xi = A) < 1/2) = 1 - g. (7) 

Then, by using the Bayes' theorem, we update agent i's belief Pr(xi = A|cr) on the basis 
of the old belief Pr(xi = A) (= 1 — Pr(xj = B)) and the perceived signal (i.e., a or /3). 
The perceived signal may be different from the received signal (i.e., a or b) because of their 
confirmation bias (Eqs. (j4])-(|5])). We assume that agents are not aware that they may be 
subject to the confirmation bias. Agents use the subjective conditional probabilities given 



by 



Pr(a 



a\s = a) = Pr(cr = /3\s = b) = 1 



(8) 



and 



Pr((T 



a\s = b) = Pr((T = = a) = 



(9) 



to perform the Bayesian update 



The posterior belief Pr(xi = A\a) is given by 




Pr(cr|xi = A) Pr(xi = A) 



Fi{a\xi = A) Pr(xi = A) + Pr(cr|xi = B) Pr(xi = B) 
E.=a,6P?(^l^) = A) Pr(a;, = A) 



, Pr(a|s) Pi{s\xi = A) Pr(xi = A) + J2s=a,b Pr(cr|s) Pr(s|x, = B) Pr(a;, = B) 




(10) 
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It should be noted that Pr(xj = B|o") = 1 — PT{xi = A\a). Then, we increment the time by 
such that each agent is updated once per unit time on average. Iterative application of 
Eq. ^ leads to 



Pr(xi = A) 



0■nc,^(l - 0)ni3i + (1 - 0)nc.iQngi 



and 



Fr(x- = = — — - 



6 



1 ^ \ '^pi~'^a.i 
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(11) 



(12) 



where n^j (n^i) is the number of signals a = a {a = (3) that agent i has perceived. The state 
of each agent i is uniquely determined by n^i — ripi, which is consistent with basic Bayesian 
theory [34, 



III. RESULTS 



A. Setup for numerical simulations 

Unless otherwise stated, we set 1/2 < 9 < 1 and assume a neutral initial condition 
Pr(xj = A) = 1/2 {I < i < N), or equivalently, Hai = ripi {1 <i < N). The agents exchange 
signals and update their beliefs, possibly under a confirmation bias. After a transient, the 
agents believe in either opinion with a strong confidence, i.e., Pr(xi = A) ^ or 1. We 
halt a run when |nQ,j — n/3j| > Anc is satisfied for all i for the first time, where ArZc is the 
threshold. 



B. Case without the confirmation bias 



We first consider the case without a confirmation bias (i.e., g = 0). We investigate the 
dynamics of the mean belief -Pa(^) = Xlili ^fi^j = A)/A^ at time t by drawing a return 
map, i.e., -Pa(^) as a function of Pa(^ — 1) (lO, |37|, l38|. The return map for = 100, 



9 = 0.99, and Anc = 500 based on 1000 runs is shown in Fig. O Because -PA(i) > -PA(i — 1) 
when 0.5 < PA(t - 1) < 1 and PAit) < Pxit - 1) when < Pxit - 1) < 0.5, the dynamics 
is in accordance with majority rule behavior. All the 1000 runs finished with an agreement 
of opinion A (i.e., Pr(xj = A) 1 for all i) or opinion B (i.e., Pr(xj = A) for all i). 
Each case occured approximately half the time. 
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C. Case with the confirmation bias 



We turned on the confirmation bias to examine the possibility that it induces disagree- 
ment among agents. At least for large q (i.e., g ~ 1), disagreement is expected to be reached, 
because the first perceived signal would determine the final belief of each agent and is equally 
likely to be a and /3 for many agents. 

In the following numerical simulations, we measured the degree of disagreement, which 
we defined as follows. We determined that agreement was reached in a run if the final signs 
of riai — ripi were the same for i = 1,. . . ,N. Otherwise, we said that disagreement was 
reached. We denoted the fraction of runs that finished with disagreement by F^. 

We set Aric = 500 and the number of runs to 1000. In Fig. [3](a) and[3]^b), is shown as 
a function of q and for = 2 and 100, respectively. First, monotonically increases with 
q and decreases with 6 for both N = 2 and = 100. It should be noted that disagreement 
occurred in at least one run in the regions right to the solid fractured lines in Fig. [3l Second, 
Fd for A^ = 2 (Fig. [3t^a)) is smaller than for A^ = 100 (Fig. Et^b)) for all the q and 9 values. 
Therefore, disagreement seems to be a likely outcome of the model for large A^, particularly 
for large q and small 6. When A^ = 100, perfect agreement, i.e., Fd = 0, is realized only for 
q close to zero. In other words, even a small degree of confirmation bias elicits disagreement 
among the agents. 



D. Probability flow analysis for N = 2 

To obtain analytical insights into the model, we performed an annealed approximation 
for A^ = 2 by averaging out fluctuations of the dynamics for different times and runs. The 
configuration of the population is specified by (mi,m2) = {riai — npi,na2 — ^^2)- The 
stochastic dynamics of the model can be mapped to a random walk on the two dimensional 
lattice; a walker is initially located at (mi, 7712) = (0, 0) and randomly hops to one of the four 
neighboring lattice points in each time step. We defined fn^rrii, 1712), /l('^i, ^2), fu{rni, 777-2), 
and /d(771i, 7712) as the probabilities that the walker located at {1711,1712) moves to {irii+l, 7712), 
(7771 — 1, 7^72), (777.1, ^2 + 1)5 and (777i, 7772 — 1), respectively. The four probabilities are given 
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by 



/fl(mi, 777,2) = < 



' (1 - g) Pr(s = a\m2) + q 
2 

Pr(s = a\m2) 



[1 — q) Pr(s = a\m2) 



K > 1), 

(mi = 0), 
(mi < -1), 



(13) 



/l(777i, m2)^ -- /ii(mi, ma), 



(14) 



fu{mi,m2) 



' (1 — g) Pr(s = a|mi) + q 



Pr(s = a|r?7i) 
2 

(1 — q) Pr(s = a|777i) 



(m2 > 1), 
(m2 = 0), 
(m2 < -1), 



(15) 



and 



where 



/i3 (7771, 7772) = - - /[/(mi, 7772), 



Pr(s = a\mj) ^ 1 - 6 + (26 - 1) Pr(xj = A|777j) 

1 - e + (2^ - 1) a + ' ^ 



(16) 
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(17) 



is the probabihty that agent j with 7iaj — nj^j — rrij imparts signal s — a. /r (mi, 7772) and 
/t/(mi, 7772) increase with mi and 7772, and fL{mi,m2) and /^(mi, 7772) decrease with 7771 and 
m2. 

In the following, we study the mean dynamics of the random walk driven by the drift 
terms. Because the transition probability of the random walk is symmetric with respect 
to the two lines 777.1 = m2 and mi — — m2, we focus on the region given by —m2 < mi < 
777.2, m2 > 0. We define 777ic and 77i2c, which are not integers in general, as the values 
satisfying fu{mic,m2) = /D(mic,m2) for ^7712 > and /^(mi, 7772c) = /L(mi,7n2c) for 
^777.1 < 0, respectively. They are given by 



7772c = -mic 



In 



{29 - 1) 



2(1 -g) 



[1-9) 



In 



9 



(18) 



1-9 
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Note that mic and exist if and only if {26 — 1) {1/2(1 — g) — (1 — 6)} ^ — 1 > 0, i.e., 

.<l-^. (19) 

First, we consider the case q < 1 — 1/26. We partition the upper quadrant of the lattice 
(given by — m2 < mi < 1712, m2 > 0) into the five regions: region 1 (0 < mi < 1712), region 2 
(mi = 0, 7712 > 0), region 3 (— m2 < mi < 0, < m2 < m2c), region 4 (mic < mi < 0, m2 > 
m2c), and region 5 (— m2 < mi < mic), as shown in Fig. IH^a). We obtain from the condition 
1/2 < ^ < 1 

fnimi, m2) - fiimi, 1712) > /[/(mi, m2) - fnimi, 1712) > (20) 



in region 1, 



in region 2, 



q 

/R(mi,m2) - /L(mi,m2) > 0, fu {thi, 7122) - fnimi, 1712) = - (21) 



fu{mi,7n2) - fD{mi,7n2) > /l("^i,"^2) - /i^("^l,"^2) > (22) 



in region 3, 

/^(mi, m2) - /L(mi, m2) > 0, /c/(mi, m2) - /d("^i, m2) > (23) 
in region 4, and 

/^(mi, m2) - /L(mi, m2) > /^(mi, m2) - /{/(mi, m2) > (24) 

in region 5. The probability flow of the walker after the annealed approximation, i.e., 
(/^(mi, m2) — /L(mi, m2), /{/(mi, m2) — f 0(^1,^2)) inferred from Eqs. (l20l) - (!2^ is shown 
schematically in Fig. ID^a). If the walker is in the second quadrant (i.e., regions 3, 4, and 
5) where the two agents disagree with each other, the random walker is likely to eventually 
escape and enter the first quadrant (i.e., region 1) where the two agents agree with each other. 
In fact. Fig. ^a), which shows the actual probability flow, indicates that the agreement 
necessarily occurs. Therefore, agreement is the expected outcome when q < 1 — 1/26. 

Second, if g > 1 — 1/26', regions 4 and 5 are absent because mic and m2c diverge. Regions 
1 and 2, in which inequalities (|20|) and (|2T|) are satisfied, respectively, are the same as those 



in the case g < 1 — 1/29. Region 3, in which inequahty ( 122|) is satisfied, is modified to 
— m2 < mi < 0. The probabihty fiows are schematically shown in Fig. HJ^b). /L(?Tii,m2) > 
/i?("^i;'^2) and fu{mi,m2) > f 0(^1,1712) are satisfied in region 3. Therefore, once the 
walker is deep in the second quadrant, it is likely to move toward mi — )■ — 00 and 1712 — ?■ 00, 
which implies that two agents finally disagree. The actual probability fiow shown in Fig. 
[5]^b) is consistent with this prediction. 

The transition line q = 1 — 1/26 is shown by the dashed line in Fig. [31(a). It accurately 
predicts the parameter region in which disagreement can occur, i.e., the region right to the 
dashed fine. 

The same transition line is also derived for the Rabin-Schrag model, which is concerned 
with a single agent subjected to a confirmation bias j3^]. In their model, the agent forms 
a belief by repetitively receiving a stochastic signal s G {a, b} from nature, according to 
Pr(s = a\x = A) = (s = b\x = B) = 9. Rabin and Schrag calculated the probability 
that the agent eventually misunderstands the state of the nature (i.e., A or B), starting 
from neutral belief. This probability is equal to zero when q < 1 — 1/29 and positive when 
q> 1 — 1/29 (see proposition 4 in [34]). Our results obtained in this section are consistent 
with theirs because disagreement in our model roughly corresponds to misunderstanding in 
the Rabin-Schlag model. 

E. Different disagreement configurations for large N 

In general, there are — 1 disagreement configurations, as distinguished by the number 
of agents that finally believe in opinion A, which ranges from 1 to A^ — 1. To distinguish 
different disagreement configurations, we examined the fraction of agents that believed in 
the minority opinion at the end of a run. We averaged this fraction over the runs ending 
with disagreement. We called this quantity the average size of the minority. 

Figures [6](a) and [6](b) show the average size of the minority for A^ = 10 and A^ = 100, 
respectively. The black regions indicate the parameter values for which the average size of 
the minority is undefined, because all the 1000 runs end with agreement. When q is small, 
the average size of the minority monotonically decreases with q and monotonically increases 
with 9 for both A^ = 10 and 100. Therefore, small q and large 9 values allow only balanced 
disagreement configurations, in which the numbers of the agents believing in the opposite 
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opinions are close to N/2. 

However, the average size of the minority increases when q is large. This is particularly 
the case for = 100 [Fig. [6]^b)]. This increase occurs for the following reason. With a 
strong confirmation bias, agents end up with an opinion consistent with a small number 
of signals perceived in the early stages, and both signals are equally likely to be observed 
in the early stages under neutral initial conditions. In the extreme case in which g = 1, 
agents reinforce the opinion that is consistent with their first perceived signal. Therefore, 
unbalanced disagreement configurations are rarely realized when q is large. 

F. Effects of the system size and initial condition 

Figures [3] and [6] suggest that the agreement is unlikely to be reached in a large popu- 
lation. To examine the effect of the population size, we defined qc as the value of q such 
that a threshold number of runs among 10^ runs end with agreement. For a given 6 value, 
we determined qc by the bisection method. The number of agreement runs may not mono- 
tonically change in q because the number of runs is finite. Therefore, the bisection method 
does not perfectly work in general. However, we corroborated that the following results were 
negligibly affected by the lack of monotonicity. 

The dependence of gc on is shown in Fig. El^a) for three threshold values. For example, 
the results for the threshold value 100 (shown by circles) indicate that at least 100 runs 
among the 10^ runs end up with disagreement when q > qc- We set 6 = 0.99 and Auc = 500. 
To explore the possibility of disagreement in large populations, we set 9 close to 1. It should 
be noted that Fig. |3] indicates that the probability of disagreement is small for a large 9 
value. In Fig. [D^a), qc quickly decreases for < 10 and gradually decreases for > 100. 
The results are consistent with those shown in Fig. |3]^b). Disagreement often occurs for 
large A^ unless q is small. Nevertheless, Fig. [Tl^a) suggests that the range of q for which 
agreement always occurs survives for diverging A^. 

In generating Fig. [Tl^a), we used an initial condition in which all the agents had a neutral 
belief (i.e., Pr{xi = A) = 0.5, 1 < i < N). To check the effect of the initial condition, we 
investigated the dependence of qc on A^ under two other initial conditions. In the bimodal 
initial condition, we initially set Uai — nj^i = 100 {I < i < N/2) and nai — npi = —100 
(A^/2 + 1 <i < N). We assumed that A^ was even for this initial condition. In the so-called 
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most unbalanced initial condition, we set Uai—npi = 100 and nai—njn = —100 (2 < i < A^). 

The numerical results for the two initial conditions are shown in Figs. Wljo) and[7|^c). 
The parameter values 6 = 0.99, An^ = 500 are the same as those used in Fig. [7](a). The 
transition point Qc decreases with more rapidly with the bimodal initial condition [Fig. 
nih)] than with the neutral initial condition [Fig. El^a)]. This result is intuitive: the bimodal 
initial condition paves the way to disagreement. In contrast, Qc under the most unbalanced 
initial condition is almost constant near 0.5 irrespective of A^. Therefore, disagreement is 
highly unlikely unless the confirmation bias is strong (i.e., q is greater than 0.5). The results 
shown in Fig. [7] suggest that the eventual behavior of the model strongly depends on the 
initial condition even after the results are averaged over runs. 



IV. DISCUSSION 

Our numerical results are summarized as follows. When the confirmation bias is absent 
(i.e., g = 0), the opinion dynamics under the Bayesian update rule leads to the complete 
agreement among agents. The behavior of the model is similar to majority rule dynamics 
(Fig. H]). When the confirmation bias is present, disagreement is a likely outcome, par- 
ticularly for a strong confirmation bias (i.e., large q). Disagreement is also more likely for 
a lower fidelity of the signal (i.e., 6 ~ 1/2) and a larger system size. The transition line 
separating the parameter region in which both agreement and disagreement can occur and 
that in which only agreement occurs is approximately given by g = 1 — 1/26 when N = 2. 
This line is identical to the one determined by Rabin and Schrag for their model for a single 



agent's decision-making 3J|. Finally, the behavior of the model strongly depends on the 
initial condition. 

35|, although Orlean's model employs 



Our model and results are different from Orlean' s 
multiple agents that perform the Bayesian updates under a confirmation bias. First, the 
belief of each agent is binary in Orlean's model, whereas our model introduces an infinite 
range of discrete beliefs, as in [34|. Second, interaction between agents is introduced dif- 
ferently in the two models. In Orlean's model, each agent refers to the global fraction of 
agents believing in one of the two opinions. In our model, agents refer to other opinions by 
peer-to-peer interaction, i.e., by receiving a binary signal that is correlated with the belief of 
the sender. Third, the stochastic dynamics of Orlean' s model is ergotic when the collective 
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opinion does not reach agreement. The collective opinion obeys a stationary distribution, 
irrespective of the initial condition. In contrast, in all our simulations, the stochastic dy- 
namics of our model was non-ergotic, such that the final configuration depended on the 
initial condition in a wide parameter region. 

In social science studies of polarization, several authors analyzed Bayesian models in 
which different agents receiving a series of common signals end up in disagreement. The 



proposed mechanisms governing disagreement include different initial be 



26 



- l30j , different update rules 



iefs or factors that 
2^], and ambiguity 



affect perception of later incoming signals 

aversion 27|, |39|. These models and ours are different in three major ways. First, a ground 
truth opinion corresponding to the state of nature is assumed in these models but not in ours. 
Second, public signals commonly received by different agents are assumed in these models but 
not in ours. Third, the agents do not have direct peer-to-peer interaction in these models 
but they do in ours. Models with interacting Bayesian agents, which show disagreement 



(reviewed in Ref. l29l), are also different from our model in the first respect. It should be 

n 

noted that Zimper and Ludwig discussed confirmation bias with their Bayesian model |27[. 
However, they derived a confirmation bias from their model, rather than assuming one, such 
that their results pertaining to confirmation bias were also distinct from ours. 

Extending our model to the case of networks is straightforward. For example, we can 
select a recipient of the signal with probability and then select the sender with equal 
probability among the neighbors of the recipient on the network. Another possible update 
rule is to select the sender first and then the recipient among the sender's neighbors. Yet 
another possibility is to select a link with equal probability and designate one of the two 
agents as sender and the other as recipient. On heterogeneous networks, the results may 
depend on the update rule because it is the case in the voter model [7I-I9I.I40I. Extension of the 
model to the case of confirmation bias heterogeneity may also be interesting. Neurological 



The 



evidence shows that different individuals have different confirmation bias strengths 
strength of the confirmation bias and the position of the node in a social network may be 
correlated and affect the dynamics. It is also straightforward to extend the model to the case 
of multiple opinion cases. These and other extensions, along with the study of analytically 
tractable models that capture the essence of the present study, are warranted for future 
work. 
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j's belief 



signal signal 
innparted by j perceived by i 



FIG. 1. Schematic of signal conversion. Signals a and a correspond to A. Signals b and /3 
correspond to B. 9 is the reliability of the signal, and q is the strength of the confirmation bias. 
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^a(^-1) 

FIG. 2. Return map of the mean belief. We set N = 100, g = 0, and 6 = 0.99. We recorded 
the values of (PaC* - 1), PaC*)) for t = 1, 1 + 1/N, 1 + 2/N, ... for 1000 runs and divided the 
recorded pairs into 21 classes. The kth class {k = 1, . . . , 21) was composed of the pairs satisfying 
{k - 1)/21 < Pxit - 1) < A;/21. We obtained the mean value (PA{t)),^ for the kth class by 
averaging PA{t) over all the pairs contained in the kth class. Finally, we plotted (Pa(*))j. against 
{k — 0.5)/21 for A; = 1, ... , 21. The diagonal PA{t) = PA{t — 1) is also shown as a guide. 
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FIG. 3. Fraction of disagreement (a) N = 2. (b) N = 100. Solid lines represent the 
boundary between i<d = and i<d > 0. The dashed line represents q = 1 — 1/29. In (a), the two 
lines almost overlap each other. The initial belief of each agent was assumed to be neutral (i.e., 
Fi{xi = A) = 0.5, 1 < i < A^). 
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(a) (b) 

TTli = —1712 ^2 m-L =7722 171-1 = ^2 = 7712 




FIG. 4. Schematic of the probabihty flow given by fR{mi,m2) — fL{mi,m2) and /[/("ii,m2) — 
/£)(mi,m2). (a) g < 1 — 1/20. (b) g > 1 — 1/20. The labels from 1 through 5 correspond to the 
five regions. 
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(a) 



(b) 




FIG. 5. Probability flow of the opinion dynamics when N = 2. Vector (/R(mi,m2) — 
/^('ni, m2), /[/(mi, m2) — fD{rni,m2)) is shown by an arrow of proportional size at each posi- 
tion of the random walker (mi, 7712). (a) q = 0.15 and 6 = 0.64, which satisfies q < 1 — 1/26. (b) 
q = 0.4 and = 0.64, which satisfies q > 1 — 1/20. The size of the vectors is manually normalized 
for clarity, independently for the two panels. 
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(a) (b) 




<7 Q 

FIG. 6. Average size of the minority, (a) N = 10. and (b) N = 100. The initial behef of each 
agent is assumed to be neutral (i.e., Pr(xj = A) = 0.5, 1 < i < N). The black region represents 
the case where all the 1000 runs end with agreement such that the average size of the minority is 
undefined. 
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(c) 

0.55 



Qc 0.5 



0.45 



10° 10^ 10^ 10^ 10"^ 

N 



FIG. 7. Threshold for the agreement-disagreement transition (gc) as a function of N under the 
(a) neutral, (b) bimodal, and (c) most unbalanced initial conditions. We set 9 = 0.99. 
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